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ABSTRACT 

We derive the luminosity function of high-redshift Lyman a emitting sources from a deep, blind, 
spectroscopic survey that utilized strong-lensing magnification by intermediate-redshift clusters of 
galaxies. We observed carefully selected regions near 9 clusters, consistent with magnification factors 
generally greater than 10 for the redshift range 4.5< z <6.7. Eleven emission-line candidates were 
located in the range 2.2< z <5.6 whose identification we justify as Lyman a, in most cases via further 
spectroscopic observations. The selection function we constructed for our survey takes into account 
our varying intrinsic Lyman a line sensitivity as a function of wavelength and sky position. By virtue 
of the strong magnification factor, we provide constraints on the Lyman a luminosity function to 
unprecedented limits of 10 40 erg s _1 , corresponding to a star-formation rate of 0.01 M Q yr~ 4 . Our 
cumulative z ~ 5 Lyman a luminosity function is consistent with a power law form, n(> L) oc L^ 1 
over 10 41 to 10 42 5 erg s _1 . When combined with the results of other surveys, limited at higher 
luminosities, our results suggest evidence for the suppression of star formation in low-mass halos, as 
predicted in popular models of galaxy formation. 

Subject headings: galaxies: formation, evolution, high-redshift, luminosity function — cosmology: 
observations — gravitational lensing 



1. INTRODUCTION 

The epoch of cosmic reionization, when the intergalac- 
tic hydrogen in the universe transitioned from neutral 
to ionized, is the current frontier of observational cos- 
mology. QSOs discovered by the Sloan Digital Sky 
Survey (SDSS) in dicate that reioniza t ion was just fin- 
ishing, at z ~ 6 llBecker et all 1200 It iDiorgovski et all 
l200lHFan et al.ll2002D . Recent results from the WMAP 
satellite suggest that signific ant reionization of th e uni- 
verse took place by z ~ 12 (jSpergel et all 12003(1 . The 
sources that reionized the universe, however, are still un- 
known: at z ~ 6 neither bright QSOs discovered by SDSS 
<Fan et al.ll200lH nor faint AGN from deep x-ray obser- 
vations l|Barger et all20 031 produced enough photons to 
reionize the universe. Other evidence from the temper- 
ature and ionization state of the intergalactic medium 
(IGM) suggests that, though QSOs dominated the meta- 
galactic ionizing backgroun d at z ~ 3, the spectrum was 
softer at reionization (e.g. iSokasian. Abel, fc HernauistJ 
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2003). Accordingly, hot stars in early star-forming sys- 
tems may have been the dominant source of reionizing 
photons. One goal of the forthcoming NASA/ESA James 
Webb Space Telescope (JWST), a 6-meter IR telescope 
scheduled for launch in 2010, is to study the formation 
of the first gene rations of galaxies and thei r contribution 
to reionization l|Mather fe Stockmanll200Ci|) . 

Early galaxies played many important roles beyond 
their involvement with reionization. The IGM was en- 
riched well above the primordial metal ab undance by 
z = 5 l|Songailal I200H iPettini et alJ 12003(1 : additional 
evidence for early metal production comes from metal- 
poor globular clusters in the Milky Way. Age estimates 
imply a formation epoch of z > 4 for curre nt cosmo- 
logical parameters llKrauss fc Chaboverl 12003). but the 
typica l metallicity of these objects is 10~ 2 of the solar 
value (Harris 1996). The stars responsible for reioniza- 
tion and early metal production may still be present in 
some form today. It is an important challenge to identify 
the transition between the very first, metal-free, stars, 
and Population II stars, because of the strong constraints 
on the metallicity of low-mass stars provided by studies 
of halo stars in the Milky Way. A complete understand- 
ing of the metallicity distribution of old Galactic stars 
will benefit from direct observation of very high redshift 
star formation in proto- galactic systems that will evolve 
into galaxies like the Milky Way. 

In advance of JWST, which will use IR capabilities to 
observe galaxies before the end of reionization in rest- 
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frame UV and optical light, current ground-based facil- 
ities have the opportunity to discover and characterize 
star-forming galaxies near the epoch of reionization with 
rest-frame UV observations. In particular, the identifi- 
cation of Lyman a emission from star-forming regions of 
early galaxies has proven to be a powerful tool for discov- 
ering z > 4 galaxies and measuring their redshifts (see 
Section 0) . The rcdshift range 5 < z < 7 is of partic- 
ular interest for two reasons. One is that the very de- 
tection of Lyman a emission may place a constraint on 
the progress o f reioniza tion llRho ads fc Malhotral 120011: 
iHaimanlhool IHu et al] l2~002at iRhoads et al] 12003ft : it 
is difficult to observe Lyman a emission from galaxies 
embedded in a neutral IGM, but the strength of the con- 
straint derived from the successful detection of L yman a 
depen ds on the assumed properties of the sources (Santos 
2003) . The second important reason to study galaxy for- 
mation during and after reionization is that an intense 
UV background and 10 4 K IGM is predicted to sup- 
press star^jojmTB^tionjn after reioniza- 
tion llBarkana fc Loeblll999t lGnedinll2000t iBenson et alJ 
12002^ 7 There is a discrepancy between some theoret- 
ical predictions of the abundance of dark matter halos 
on dwarf-galaxy mass scales and the number of observe d 
dwarf satellites in the Local Group l|Moore et alj ri999) . 
Reionization may sterilize many dwarf galaxy-scale halos 
to star formation, so that the luminous satellites of the 
Milky Way are dwarf galaxies formed before reionization, 
and the re maining satelli te halo s are empty of stars and 
thus dark l|Benson et alJl2002ri ISomervillell200l . 

This paper presents the results of a spectroscopic 
Lyman a emission-line survey for galaxies at 2.2 < 
z < 6.7 that utilizes the strong- lensing properties of 
intermediate-redshift clusters to magnify the surveyed re- 
gions. In Section|21we review the use of Lyman a surveys 
as probes of early star formation. Section |31 motivates 
the importance of surveys for low-luminosity galaxies. 
Section 2| describes the advantages of a survey utilizing 
strong lensing and details our survey strategy, targets, 
observations, and data reduction. In Section our Ly- 
man a emission-line detections are presented. We com- 
pute our survey volume and source number density in 
Sectional Section compares the results of our survey 
to other surveys and theoretical models. In Section|Hlwe 
summarize. 

Throughout this paper we use a ACDM cosmologi- 
cal model with (fi m , fi A , ^b, cr 8 )=(0.3, 0.7, 0.043, 0.9) and 
h = Jf /(100 km s~ 1 Mpc _1 ) = 0. 7; these values are co n- 
sistent with the values derived in lSpergel et al.1 <)2003ft . 

2. LYMAN a SURVEYS 

The primary appeal of Lyman a emission as a sign- 
post to high-redshift galaxy formation is that it traces 
star formation at a wavelength that conveniently red- 
shifts into the visible and near-IR, where sensitive, high- 
angular resolution observations are currently most prac- 
tical. The Lyman a emission line may be quite strong, 
but its luminosity is quite sensitive to physical details of 
the nature and geometry of the star-forming regions. Be- 
cause Lyman a emission traces hydrogen recombinations, 
it is intimately related to the production of ionizing pho- 
tons by the stars present. Both the initial mass function 
(IMF) and metallicity of the stars affect their production 
rate of ionizing photons. However, if the IMF and metal- 



licity are constrained to reasonable ranges, the ionizing 
photon production rate can be reliably connected to the 
star-formation rate. 

The major complication for the interpretation of Ly- 
man a line strengths is the effect of the nebula around the 
star-forming region. Hydrogen at low density does not 
recombine quickly, so, e.g., ionizing photons that escape 
into the IGM are "lost" for the purposes of producing a 
Lyman a emission line. Even after a hydrogen recom- 
bination produces a Lyman a photon, whi ch happens 
for a bout two thirds of the recombinations ijOsterbrockl 
11989ft . there are many ways in which the photon may 
be destroyed prior to escape. The resonant nature of 
the Lyman a transition results in a very short mean free 
path, even in a mostly ionized nebula. Consequently, if 
dust is mixed in with the gas, then the chance of ab- 
sorption by a dust grain may be higher for a Lyman a 
photon than a non-reson antly scattered photon at the 
same wavelength (but see lNeufeldlll99lL for an alterna- 
tive). The dust content of very high-redshift galaxies is 
still relatively unconstrained, and will likely remain so at 
least until JWST. 

On the positive side, Lyman a is the intrinsically 
strongest recombination line from an H II region. An- 
other meritorious aspect of Lyman a is that its emission 
strength does not strongly depend on the metallicity of 
the ionized gas (the only effect is from the temperature 
of the photoionized gas, which depends on metallicity); 
consequently, it can be used as a tracer for truly primor- 
dial star formation, where dust extinction is also believed 
n ot to be a problem. 

iPartridee fc Peebles! l)1967ft introduced a model of 
galaxy formation "to assess the general possibility of 
observing distant, newly formed, galaxies." In their 
model they estimated that 6-7% of the luminosity of 
early galaxies would be emitted in the Lyman a emission 
line: they predicted line luminosities of 2 x 10 45 erg s _1 
over galaxy f ormation time-scales of 3 x 10 7 yr. The 
predictions of IPartridee fc Peebles! l)1967ft and others led 
to many observational surveys for Ly man a emission 
from hi gh-redshi ft galaxies, reviewed bv iPritched (1994, 
§4.5). Pritchet summarized the status of searches at 
that time with "no emission line primeval galaxies have 
been found" despite 16 cited surveys covering various 
redshifts ranges from z = 2 to z = 5. In striking contrast 
to these pioneering explorations, many high-redshift Ly- 
man a emission-line galaxies h ave been discover e d and 
confirme d in the past 8 years. IStern fc Soinradl l)1999L 
§4) and iTaniguchi et al.1 l)2003ft have reviewed recent 
progress, and we provide a brief summary here. 

The search technique that has dominated recent suc- 
cess in the discovery of large numbers of Lyman a 
emission-line galaxies at z > 4 is narrow-band pho- 
tometry. This approach uses a narrow (~ 100 A) 
filter chosen to lie in a spectral region of low sky 
background; such surveys cover relatively large areas 
of sky with sensitivity to Lyman a emission over a 
small window in rcdshift, Az ~ 0.1. Many groups 
have now performed su ccessful blind narrow-band sur- 
veys for z > 4 gala x ies: IHu. Cowie. fc McMahonl 1119981) 
and IRhoads et al J KOOffl at z = 4.5: lOuchiXaTl 
ll2003ft at z = 4.9: IHu McMahon. fc Cowie] lll999ft . 
IRhoads fc Malhotral f2001ft. a nd IRhoads et all (12003ft at 
z = 5.7: and IHu et alJ l|2002aft and lKodaira et alj '(2003) 
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at z = 6.5. 

Spectroscopic surveys provide a complementary tech- 
nique to the narrow-band method. For equal observ- 
ing time at one position, spectroscopic searches at op- 
tical wavelengths can cover a large range in redshift, 
Az ~ 3, to better line flux sensitivity than a correspond- 
ing narrow-band survey. However, the area surveyed by 
a long slit is typically ~ 5 x 10~ 2 arcmin 2 , in contrast 
to 20-2000 arcmin 2 for an imaging camera. Deep long 
slit observations of other targets have been searched for 
serendipitous detection of high-rcdshift galaxies. These 
techniques disco vered the first confirmed z > 5 galaxy 
ijDev et alll998l) and subs equen t ly turned up a few more 
z > 4 sources llHu et allll99l iWevmann et alJll998jh 
Serendipitous surveys will likely continue to play a role 
in discovering high-redshift galaxies, since the deepest 
spectra on the largest telescopes are likely to be pointed 
observations rather than devoted emission-line surveys. 

Under sky-limited conditions, it is simple to show 
that the signal-to-noise ratio 1Z that is reached on pure 
emission-line sources of flux F in a survey over sky area A 
and wavelength range A, with area and wavelength cov- 
erage AA and AA, seeing disk S A and spectral resolution 
SX, is given by 

n= Ft ^ F^/Tj{AK) 

y/SSA SXt I q SA SX ' 

y AA AA 

where t and T are the lengths of individual exposures, 
and of the whole survey, respectively, and S is the sky 
surface brightness per unit wavelength. The tradeoff is 
thus between the number of seeing elements that can 
be observed simultaneously (0(N 2 ) in case of an imag- 
ing survey, O(N) for a long-slit spectral survey) vs. the 
number of spectral elements (1 for a narrow-band survey, 
0(N) for a spectral survey). As each is limited by detec- 
tor sizes, both modes can in principle offer comparable 
survey speed to a given flux limit (the argument applies 
also to surveys with integral field unit spectrographs). 
Spectral surveys are advantageous for targeting small re- 
gions of the sky and covering a large wavelength range, 
or for simply going as deep as possible in a given amount 
of time; narrow-band surveys are optimal in the opposite 
regime of wide area, small wavelength range. A further 
advantage of high-spectral resolution surveys is that the 
night sky emission lines can be maximally avoided, which 
is not possible with the typical bandwidths employed in 
narrow-band surveys. 

A few other hybrid techniques combine aspects of 
the narro w-band imaging and long-slit spectroscopy ap- 
proaches. iMaier et al.1 1)20031) used an imaging Fabry- 
Perot interferometer to take sequences of narrow band 
images within the night-sky windows corresponding to 
z = 4.8,5.7,6.5, and have c onfirmed d is coveri es at 

z = 4.8 and 5. 7. Recently iLillv et alJ l|2003D and 

iMartin fc Sawickil l)2003|) used a "slit-let" slit mask with 
a narrow-band filter to do spectroscopic surveys over 
a relatively large areas and narrow redshift windows. 
There is also an ongoing HST program to use slitless 
spectroscopy with the ACS camera grism to discover 
high-redshift Lyman a emission (J. E. Rhoads, priv. 
comm.). 

3. SEARCHES FOR DISTANT LOW-LUMINOSITY 
GALAXIES 



In this section we present the motivation for conduct- 
ing a survey devoted specifically to low-luminosity z ~ 5 
Lyman a emitting galaxies (< 10 42 erg s -1 in the Ly- 
man a line), an unexplored region of survey parameter 
space. 

At z ~ 5, t he luminous Lyman a galax ies (0 
and QSOs (e.g., iFan et all lillM EM 120011) discov- 
ered so far represent the rarest and most spectacular 
tail of the range of structure formation scales (e.g., 
iBarkana fc Loebll2003l suggested that the high-z SDSS 
QSOs reside in 10 12 M Q virialized halos). They almost 
certainly evolve into the rarest and most massive envi- 
ronments in the local universe: rich clusters of galaxies. 
In striking contrast, the characteristic mass of virialized 
halos at z ~ 5 is only 10 9 M ; if such an object steadily 
converts its ~ 10 s M Q of baryons into stars, it will have 
a star-formation rate of only ~ 0.1 M yr _1 . 

These objects, which we will refer to as low-mass halos, 
would not be detectable in any of the surveys cited above, 
but our understanding of galaxy formation depends cru- 
cially on constraining their properties for three impor- 
tant reasons: First, they represent the most common 
environment by mass of virialized halos (the peak of the 
mass-weighted mass function is always near the charac- 
teristic mass); consequently, if low-mass halos form stars 
efficiently, they could dominate the star-formation rate 
at high redshift. Second, they are the progenitors of 
common galaxies in poor environments, like the Milky 
Way, under the current paradigm of structure forma- 
tion. The detection of low-mass sources is a direct test of 
the "bottom-up" description of galaxy assembly. Third, 
low-mass objects have a unique link to the IGM: ki- 
netic energy injected into the IGM by photoionization 
at z > 6 is expected to raise the cosmic Jeans mass 
and inhibit gas cooling; this effect has no consequence 
for the brightest sources residing in deep potential wells, 
but may h eavily suppress star formation in 10 9 M^ ob- 
jects (e.g., |BarSMTLoe3[l99l [Gnedii3|200(J. This 
Jeans-mass effect has been cited as the solution for the 
Cold Dark Matter "crisis" of over-predicting the number 
of Milky Way satellite galaxies compa red to observation 
ijBenson et al.ll2002bl ISomervillell2002|) . Low-mass halos 
also place relevant constraints on the energy scale of dark 
matter in Warm Dark Matter models (Z. Haiman, priv. 
comm.). 

A practical and strategic advantage in characterizing 
the luminosity function at low luminosities is to deter- 
mine the optimum survey depth for future surveys that 
aim to discover large numbers of 4 < z < 7 galaxies. 
If the luminosity function were very steep, then deep 
surveys such as ours would be more efficient than shal- 
lower, wider field surveys. Theoretical prejudice suggests 
that the luminosity function should have a steep effective 
slope in the region associated with the exponential cut- 
off of the number density of underlying halos, assuming 
some sort of mass-to-luminosity correspondence. Previ- 
ous non-detectio ns combined w ith recent successes seem 
to bear this out l|Pritchetlll994l and see below). That is, 
current surveys may be approaching the characteristic lu- 
minosity; however, the characteristic luminosity and lu- 
minosity function shape have yet to be well constrained. 
In the local universe, luminosity functions based on star- 
formation rate estimators such as Ha luminosity roughly 
follow the Schechter function l)Schechterl 119761) form of 
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power-law behavior at low luminosity, reflecting the un- 
derlying power-law of the mass function, albeit with pos- 
sibly a different slope (e.g.. lGallego et aLlll995|) . 

The predicted suppression of star formation in halos 
with small potential wells suggests that the Lyman a lu- 
minosity function at z <~ 5 may have a modified shape. 
The reason is that the characteristic mass scale where 
the exponential tail and power law regions of the mass 
scale meet, ~ 10 9 M Q , corresponds quite closely to 
the mass scale where a number of physical mechanisms 
may suppress star formation. We described the effect 
of a hot IGM above. In addition, energy injected in 
the ISM of star-forming galaxies by stellar winds and 
supernovae (called "negative feedback," or just "feed- 
back") is predicted to heavily suppress star formation 
in halos with circula r velocities below about 100 km s _1 
ijDekel fc Silld 11986^ . corresponding to a mass scale at 
z ~ 5 of ~ 10 11 Mq. A curren t implementation of feed- 
back by iBenson et alJ l)2002al) shows the importance of 
feedback on the high-redshift star-formation rate. This 
effect complements the inhibiting effects of a hot IGM 
on star formation in low-mass halos. A third mechanism 
that may reduce star formation preferentially in low-mass 
halos is the effect of large-scale winds blown by star- 
formin g g alaxies. IScann apieco. Ferrar a. fc Broadhurstl 
(2000) and IScannapieco fc B roadhurstl l|2001f) computed 
the influence of winds blown out through the IGM by the 
first galaxies to form. They concluded that these winds 
effectively sweep gas out of nearby halos in the process 
of collapsing, meaning that even though the dark matter 
continues its collapse to virialization, there is little corre- 
sponding star formation because of the lack of baryons. 
In their model winds influence the entire star-formation 
history of the universe, but at z ~ 5 particularly suppress 
star-formation in halos smaller than ~ 10 10 M Q . 

Figure^is a schematic illustration of the possible effect 
of the suppression of star formation in low-mass halos. 
The solid curve is the mass function of halos at z = 5 
l|Sheth fc Tormenl2 002). and we have converted the mass 
scale on the top axis into a star-formation rate on the bot- 
tom axis using a simple prescription (see Section f7.2|) . In 
this simple model, the Lyman a line luminosity func- 
tion, if interpreted as directly proportional to that of the 
SFR, would have the shape of the solid curve. We have 
introduced the suppression of star formation in low mass 
halos by using the prescription of iGnedml <)2000f) to ef- 
ficiently filter out gas from halos below a critical mass 
scale Mp, the filtering mass l|Gnedinl l20001. Each of the 
broken curves is the shape of the Lyman a line luminosity 
function we expect (again assuming it scales with SFR) 
after applying filtering on a different mass scale. Ulti- 
mately the filtering mass, and thus the physics described 
in the previous paragraph, may be constrained directly 
by a measurement of the shape of the star-formation rate 
function, along with estimates of the corresponding halo 
masses (see Section l7"2*)) . 

Only recently have deep observations with large tele- 
scopes even detected z ~ 5 Lyman a emitting galaxies; 
the detailed form of the luminosity function will not be 
well-constrained in the immediate future. However, our 
low-luminosity Lyman a survey, and others like it, in con- 
cert with surveys at higher luminosities, may constrain 
or detect the break in the luminosity function associated 
with the characteristic halo mass and where star forma- 



tion is suppressed. 

4. OBSERVATIONS 
4.1. Survey Stategy 

The goal of our survey is to extend the census of Ly- 
man a sources at z ~ 5 to the faintest luminosities pos- 
sible with existing observational facilities. We achieved 
this through deep spectroscopic exposures on areas of sky 
strongly magnified by gravitational lensing. 

Strong gravitational lensing by rich clusters of galaxies 
at z ~ 0.2 is an invaluable resource to a survey covering 
a very small area of sky to great depth, as such clus- 
ters magnify background high-redshift sources by greater 
than a factor of 10 over regions of roughly 0.1 arcmin 2 in 
the image plane. This advantage comes at a price: lens- 
ing increases the apparent area of a background source at 
fixed surface brightness, so that the pointing is deeper by 
a factor of the magnification, Ai, but covers only 1/A4 
of the area of an equivalent unlensed pointing (assum- 
ing unresolved sources). This is a superior strategy for a 
deep, narrow survey. Achieving the same depth and area 
in an unlensed survey would require a factor of M. more 
observing time, a huge difference for M. ~ 10. 

Only a small area of the whole sky, ~ 100 arcmin 2 , 
is magnified by a factor of 10 or more by clusters; this 
sets a fundamental limit to the areal coverage of any sur- 
vey utilizing strong lensing by clusters. For a large sur- 
vey that detected many sources, statistical information 
about the distribution of magnification over the survey 
area might be sufficient to construct a Lyman a emitter 
luminosity function. However, since we planned to sur- 
vey a small area and detect only a few sources, we chose 
to observe clusters with HST imaging and spectroscopic 
rcdshifts f or many of the iden tified arcs and multiple im- 
ages (e.s.. iKneib et al1 fl996). These are currently avail- 
able for only a small fraction of strong- lensing clusters. 
The positions and spectroscopic redshifts of the arcs con- 
strain the distribution of mass in the cluster which can, 
in turn, be used to predict a magnification map of the 
cluster for high-redshift sources. Since lensing depends 
on the angular-diameter distance between the source and 
the lens, and that distance changes slowly with redshift 
at z ~ 5, the sky area with large magnification is fairly 
independent of the source redshift for z > 3 

Given that high-quality lens models are only available 
for about a dozen clusters observable from Hawaii, the 
total area available to us for a survey is currently quite 
small. To make the most of this limited resource, we 
conducted a spectroscopic survey. The advantage of a 
spectroscopic survey is that Lyman a emission can be 
simultaneously surveyed for over a redshift range 2.2< 
z <6.7. The primary draw-back of a spectroscopic survey 
was that, with the instrumentation available, the most 
efficient technique was slit spectroscopy. The geometry 
of a long slit is not well matched to the lensed region of 
sky (see Fig. |5J), thus some of the slit area covers area 
outside the cluster that is not strongly magnified. 

A long-slit survey does have other advantages. We can 
expect many emission-line sources in our survey other 
than Lyman a; in particular, optical lines associated 
with strong star formation, such as [O II] 3726, 3729 A; 
H^ 4861 A; [O III] 4959, 5007 A; and Ha 6563 A. Low- 
resolution spectroscopy with large wavelength coverage 
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allows the rejection of many potential low-redshift con- 
taminants through the identification of other emission 
lines. However, the [O II] doublet can be difficult to re- 
solve at low dispersion and, if redshifted to z ~ 1, there 
are often no other strong emission lines present in the op- 
tical spectrum. Thus final identification of an emission 
line as Lyman a may require follow-up spectroscopy at 
intermediate resolution. 

Redshift identification is aided by two additional fac- 
tors. Deep optical broadband imaging, available in at 
least one band for well-studied clusters, can be used, as in 
narrow-band searches, as a rejection filter: if a putative 
Lyman a system shows much observable flux shortward 
of Lyman a, then it is not likely a correct line identifica- 
tion, because the intrinsic UV spectrum combined with 
IGM absorption create a s trong decrement across t he Ly- 
man a emission line Ce.g.. lSongaila fc Cowidl20?)l . The 
second tool is available when two or more images (due to 
strong lensing) of the same high-redshift source are dis- 
covered. In this case the lensing model itself may place 
a reasonably strong constraint on the redshift of the sys- 
tem based on the o bserved image positions and flux ratios 
lIKneib et alJll99l . 

The deepest survey for a given observing time would be 
to devote all of the time to a single slit position. However, 
we expect Lyman a sources to be clustered, resulting in a 
non-Poisson distribution. To estimate an accurate lumi- 
nosity function we surveyed several independent volumes 
(via surveying behind several lensing clusters) to amelio- 
rate cosmic variance and recover the Poisson noise limit 
(see Section IT2|) . 

4.2. Survey Parameters 

Tablcnsummarizes the 9 lensing clusters of our survey. 
We have constructed a detailed mass model for each, 
based on HST imaging and lensed arc redshifts. 

Details of our spectroscopic observations are listed in 
Table [21 Clusters that were observed at multiple po- 
sition angles are designated further by an identifica- 
tion number. We used the dou ble-beam Low R esolution 
Imaging Spectrograph (LRIS, IQke et aDll995|) in long- 
slit mode on the Keck I 10-meter telescope at Mauna 
Kea to perform our survey. For the 2000 March ob- 
servations we used a slit 0.7" wide and a spectroscopic 
range of AA 6800-9500 A (corresponding to Lyman a 
with 4.6< z <6.8), using a 600-line grating blazed at 
A7500 A which gave a resolution of ~3.0 A. In 2001 April 
we switched to a 1.0"-wide slit and used a 600-line grat- 
ing blazed at 1 /jm over the same wavelength range as 
above, at a resolution of ~4.0 A. Starting in 2001 April 
we also began using a 300 line grism blazed at 5000 A and 
a dichroic at 6800 A to simultaneously take spectra on 
the blue arm of the spectrograph, over AA 4000-6700 A 
(corresponding to Lyman a from 2.2< z <4.5) at 3.5-4 A 
resolution. 

The length of the spectrographic slit was 175". We 
mapped an area on the sky via offsetting the telescope 
perpendicular to the long axis of the slit by a distance 
equal to the slit width. At each slit position we made two 
exposures of 1000 sec to facilitate cosmic ray rejection. 
Each map comprised 5-10 adjacent slit positions at the 
same position angle, giving contiguous survey areas on 
the sky of 875-1750 arcsec 2 . 

The pointing of the slit on the sky was verified by reg- 



istration of images from the LRIS slit-viewing guide cam- 
era to the HST images of the cluster (because the clus- 
ters are rich in bright galaxies, there were always many 
sources in the slit- viewing guide camera images). Our se- 
quence of slit offsets typically agreed with a regular spac- 
ing of the slit width to a precision of 0.1" (10% of the 1" 
slit). More importantly, registration of the slit position 
on the HST image enabled us to look for a broadband 
counterpart at the location of emission lines detected in 
our spectra. 

The areas mapped by the procedure above were chosen 
to take advantage of the strong magnification of back- 
ground sources provided by the foreground cluster. The 
magnification map of the sky around the cluster is con- 
strained by the distribution of visible light in the clus- 
ter and the measured velocity dispersion of some clus- 
ter members, but is crucially verified and refined by in- 
cluding information from the locations and redshifts of 
strongly lensed sources. These background galaxies, gen- 
erally brighter and at lower redshift than the z > 4.5 
sources we searche d for, have been th e target of previous 
observations fe.g.. lKneib et al.l fl9961 . 

We used up-to-date cluster mass models to gener- 
ate redshift-dependent maps of the magnification toward 
4 < z < 7 sourc es using the LENST00L software developed 
bv lKneibl ljl993j) . The angular diameter distance between 
the cluster and those redshifts depends only weakly on 
redshift, so we were able to choose areas on the sky with 
high magnification over our entire redshift range of in- 
terest. 

The geometry of the magnification map is generally 
characterized by two concentric ring-like curves of for- 
mally infinite magnification, called the inner and outer 
critical lines. These are related to the location in the 
image plane of the caustic of an elliptical potential 
(Blan dford fc Naravanl Il992(l . but modified by the de- 
viations of the projected lensing potential from an el- 
lipse. The areas of highest magnification are found next 
to the critical curves, so our survey maps generally fol- 
low the outer critical line. The outer critical line is 
more amenable to long-slit mapping because of its greater 
length on the sky and its less curved form. However, in 
one case (Abell 1689) we mapped sky near the inner crit- 
ical line as well. Lensed sources close to the critical line 
are often multiply imaged, forming a merging pair on ei- 
ther side of the critical line. We considered this when we 
chose our map locations, but the irregular shape of the 
critical line, compared to the straight shape of our slit, 
limited the extent to which we could map exclusively one 
side of the critical line. 

Figure [21 summarizes the adopted strategy for each 
cluster in the context of the location of the critical line 
for a lensed source at z — 5 (dotted lines). In the 
most massive clusters with the best mass models, such 
as Abell 1689 and Abell 2218, we made multiple maps 
(see also Table [2Jl . In these cases each survey region is 
labeled by a number corresponding to the observations 
listed in Table H 

The total area on the sky covered by our survey was 
4.2 arcmin 2 . The effective areal coverage of the survey is 
smaller due to lensing by a spatially variable magnifica- 
tion factor (see Section fo.l.lfl . 

4.3. Candidate Selection and Catalog 
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The 2-D spectroscopic data were reduced using stan- 
dard techniques in the NOAO/IRAF software environ- 
ment 4 . Cosmic rays were rejected from each pair of im- 
ages at a given loc ation with the L. A. COSMIC routine 
ijvan Dokkuml l200ll) and sky emission was removed by 
subtracting block-filtered data. Sky subtraction was not 
photon-limited on the strong night sky lines due to the 
presence of fringing features; we account for this when 
determining our sensitivity in Section [f)l We calibrated 
our absolute ef ficiency with observations of spectropho- 
tometry stars ijMassev fc G ronwalll ll990t) . 

The sky-subtracted 2-D spectral images were indepen- 
dently inspected by two of the authors (RSE and JR) 
and a catalog of 46 candidate Lyman a emission lines 
was compiled. Astrometric positions were determined 
for each and the HST images inspected for sources at 
the relevant location. In some cases, candidates were 
located beyond the boundary of the HST images and 
ground-based images were used. 

Candidate Lyman a emission lines were characterized 
on the basis of several criteria. First, the full spectrum 
(generally 4000-9500 A) was closely examined for other 
emission lines. On a second pass, candidates were ranked 
as marginal or promising depending on their apparent 
strength and spatial extent. Out of the initial list of 46, 
7 sources are confirmed Lyman a lines at 2.8 < z < 5.7, 
and 4 sources are promising candidates that we identify 
as likely to be Lyman a lines. Those 11 sources are 
listed in Table with observed line fluxes for detections 
at z > 4.5 (on the red arm of the spectrograph). 

4.4. Intermediate-resolution Spectroscopy 

The [O II] 3727 A doublet has a rest-frame separation 
of 3 A. Thus at z ~ 1, when the doublet is redshifted into 
our most important spectral range, the observed doublet 
separation is ~ 6 A. This is close to our LRIS spec- 
tral resolution of 4 A, so to determine whether any of 
our Lyman a line candidates were unresolved [O II] dou- 
blets, we followed up 15 candi dates with the Eche llette 
Spectrograph and Imager (ESI. iSheinis et alJl2002ft . We 
took spectra using the echelle mode and a 0.75 " slit, 
which delivered a spectral resolution of R — 6000. The 
exposure times varied depending on candidate strength. 
ESI spectroscopy confirmed three candidates as Lyman a 
emission lines (see Table OJ, at z = 3.27, z = 3.62 , and 
the z = 5.57 galaxy presented in lEllis etal] 1)200 If) . Ad- 
ditionally, several of the original candidates turned out 
to be [O II] 3727 A at z ~ 1. 

5. DETECTIONS 
5.1. z > 4.5 

We detected three convincing Lyman a sources with 
z > 4.5 (see Table. EJ; each was detected in photometric 
conditions. Figure |3| shows their two-dimensional spec- 
tra, and Fig. ^ shows a magnified view of the Lyman a 
emission line as well as HST images of the source loca- 
tions. 

We discovered a source tow ard Abell 2218 at z = 5.6 
that we discussed in detail in lEllis et al.1 l)2001|) . That 

4 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



source is strongly magnified (a factor of 33) and multiply 
imaged: we used the HST data to locate a second image 
outside of our survey region. The rcdshift identification 
was confirmed by an intermediate-resolution spectrum of 
both images, showing them to be identical, with P-Cygni 
line profiles characteristic of Lyman a, and certainly not 
the [O II] 3727 A doublet. Our lensing model additionally 
constrained the redshift of the source to be consistent 
only with the identification of the line as Lyman a. The 
unlensed luminosity in the Lyman a line is (7.8 ± 0.8) x 
10 41 erg s _1 . 

We blindly recovered a z = 4.89 galaxy in 
Abell 1689 that was di scove red originally by 
iFrve. Broadhurst. fe Bemtezl j2002ft. This object is 
multiply imaged, and the iFrve et al.l l)2002D spectrum 
shows a strong break across the line and metal ab- 
sorption lines in the continuum redward of Lyman a, 
confirming the redshift. We estimate the magnification 
of this source at a factor of 7.2, within th e range of 
3-14 suggested by the coarser modeling of IFrve et all 
( 2002) . The unlensed luminosity in the Lyman a line is 
(7.4 ± 0.7) x 10 42 erg s" 1 . 

In the field of C11358, we discovered a source at z = 
4.92, the sam e redshift as t hat o f the strongly lensed arc 
discovered bv lFranx et al.l 1^997). We believe our source 
is likely associated, as an additional component, with 
that responsible for the giant arc. The magnification is 
x 10, giving an unlensed luminosity in the Lyman a line 
of (2.5 ± 0.3) x 10 42 ergs" 1 . 

In addition to these confirmed sources, we discovered 
two more sources that we consider likely to be Lyman a 
emission lines at z > 4.5 (see Table. |3J). We discovered 
a source in the field of Abell 773 that we consider to be 
a good candidate for Lyman a emission at z = 4.74. If 
this is the correct identification, the unlensed luminosity 
in the Lyman a line is (2.8 ±0.6) x 10 41 erg s _1 , using a 
magnification factor of 9.5. The other likely source was 
discovered in non-photometric observations; though we 
lack absolute flux calibration for this source, we applied 
a systematic photometric correction that we consider un- 
certain up to a factor of approximately 2. The likely 
z = 4.77 source is in the field of Abell 963, magnified 
by 2.2 times. The source would then have an unlensed 
Lyman a line luminosity of (1.4 ± 0.2) x 10 42 erg s _1 
(statistical error only). Two-dimensional spectra of the 
emission lines and images of the likely candidates are 
provided in Fig. 

5.2. z < 4.5 

We detected with certainty four Lyman a emission-line 
sources at z < 4.5. One was a blind recov ery of a z = 
2.80 s ource behind Abell 370 discovered bv llvison et alJ 
(2002), and a second is another galaxy also at z = 2.80 
in the same field. Two are new detections, at z = 3.27 
and z = 3.62, both behind Abell 963. Two more good 
candidate Lyman a emission-line source identifications 
are pending. In addition two galaxies behind Abell 2218, 
both at z ~ 2.5, were detected in Lyman a absorption. 
We measured secure redshifts for 51 other sources using 
other emission lines, primarily [O II], [O III], and Ha. 
These data are useful for further constraining the cluster 
mass models, and will be presented in a separate paper 
(Richard et al., in prep.). 
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6. SURVEY ANALYSIS 

In this section we compute the number density of Ly- 
man a emission-line galaxies in our survey as a function 
of Lyman a line luminosity, L. To accomplish this we 
first determine the effective volume of the survey as a 
function of the luminosity and redshift of a source. 

6.1. Survey Volume 

A location in our survey volume is characterized both 
by location on the sky, fl, and a redshift, z. The differ- 
ential volume element located at position (Sl,z) in our 
survey is 



dV c (U,z) 



1 



M(n,z) 



d/ c (z) ( 
dz 



[D 2 c (z)dn]. (2) 



The first factor corrects for the lensing effect, which de- 
creases the area surveyed. The second factor is the co- 
moving length of the volume element along the line of 
sight, with 



dl c (z) 



Ho [^m(l 



i3/2 



(3) 



(we have assumed a flat universe). The third factor is the 
comoving transverse area of the volume element, with 



D c {z) 



dl c (z') 
dz' 



dz' 



(4) 



(subscript "c" denotes that the quantity is measured in 
comoving coordinates, which we use throughout). 

Every volume element of our survey is characterized 
by a limiting Lyman a line luminosity, Ln m (£l, z); a 
source with Lyman a line luminosity L will be detected 
in our survey provided it resides in a volume element with 
Lii m (Jl, z) < L. The limiting luminosity of a volume ele- 
ment depends on the magnification (due to lensing by the 
foreground cluster), M(£l, z), the limiting observed Ly- 
man a line flux fu m (z), and a slit transmission function, 

T(n), 



L\ im (tl, Z) = 



4n(l + z) 2 D 2 c (z) h m (z) 



r(n) 



M(n,z) 



(5) 



6.1.1. Magnification, M(fl, z) 



The magnification due to lensing by a given cluster 
is a function of both position and redshift. The refer- 
ences for the cluster mass models are given in Table 
These model s were run through the LENSTDOL software 
l)Kneiblll993j) to generate the magnification as a function 
of redshift at every position in the survey. In practice 
the area of the survey was divided into parcels of sky 
of length 0.8" and width equal to the slit width, and the 
magnification was calculated at the center of each parcel. 
The magnification at each position was sampled at nine 
rcdshifts, and the magnification at a particular redshift 
found by interpolation. 

Figure shows the magnification as a function of po- 
sition along a slit observed in Abell 2218, for two dif- 
ferent redshifts. The magnification at a given position 
is a weak function of redshift for magnification values 
less than ~ 100 (94% of the survey area), because the 
angular diameter distance between the cluster and the 
source changes by less than 25% over the source redshift 



range 4.5 < z < 6.7. Very near the critical lines mag- 
nification is a stronger function of redshift. Our survey 
maps sky by observing adjacent slit positions, so errors 
associated with interpolating the highest magnification 
A4(Q,z) > 100 points should not be important in our 
estimate of the survey volume. 

Figure0is a cumulative histogram of the magnification 
factor over the survey; the area surveyed as a function 
of magnification is very weak function of redshift, even 
at the highest magnifications. About half of the area we 
surveyed is magnified by at least a factor of 10, with the 
lower magnification values coming from area at the ends 
of the slits, because most strong-lensing clusters subtend 
a size less than the slit length (175") on the sky (see 
Fig. HI- 

6.1.2. Limiting Lyman a line flux, fn m (z) 

We define our limiting Lyman a line flux as the sig- 
nal in an aperture of 1.3" by 7.7 A that exceeds 5 times 
the root-mean-square fluctuations (noise) in apertures of 
that size, i.e., a 5-tr limit. The spatial dimension of the 
aperture was chosen to be roughly matched to the see- 
ing, and the spectral dimension was chosen match the 
expected line-width of high-rcdshift Lyman a emission 
from galaxies, ~ 300 km s -1 . If a source is larger than 
our aperture, which is especially possible along the spa- 
tial direction if the source is strongly lensed, then we will 
not be as sensitive to that source as we estimate. 

We assumed that the sky noise was constant over the 
length of the slit at fixed wavelength. This allowed us to 
include the non-Poisson contribution to the noise level 
from fringing features, which dominated the noise on 
strong sky lines. All three of our confirmed z > 4.5 detec- 
tions were more than 5-<r detections, but we found that 
visual inspection generated candidates (some of which 
were subsequently confirmed as bona fide emission lines) 
with fluxes below the 5-tr limit; in particular, one of the 
likely candidates at z > 4.5 is just at the 5-cr limit. Thus 
a 5-cr limit should be appropriate for the calculation of 
our survey volume. 

The limiting line flux varies as a function of wavelength 
due to the wavelength dependence of the atmospheric ab- 
sorption and the sensitivities of the telescope and instru- 
ment, but the largest dependence is due to the strong 
variation in atmospheric emission from OH airglow lines 
(except at A > 9300 A, where the sharp drop in instru- 
mental sensitivity dominates). Figure[H]shows fum(z) for 
the slit pointing illustrated in Fig. [U a 2000 sec observa- 
tion under photometric conditions. We compute fu m (z) 
by simple conversion of the observed wavelength into the 
corresponding redshift for Lyman a to be observed at 
that wavelength. 

Approximately half of our survey data were taken in 
non-photometric conditions. We account for this by di- 
viding the limiting line flux measured from the observa- 
tions by our best estimate of the sky transparency dur- 
ing the exposure. During some exposures we have se- 
quences of guide-camera observations that were used to 
measure relative transparency between observations, and 
in some cases absolute transparency when photometric 
guide-camera images were available. In other cases we 
rely on observation log notes based on the count rate of 
the guide star as reported by the telescope operator. 
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6.1.3. Slit transmission, T(f2) 

The slit widths used in our survey, originally 0.7" and 
later 1", are comparable in size to the seeing disk. Con- 
sequently the transverse distance of a source from the 
center-line of the slit has a small impact on the source's 
observability: objects in the center of the slit are easier 
to detect than objects at the slit edge. Since the absolute 
calibration was performed with respect to standard stars 
in the center of the slit, we compute the fraction of light 
transmitted through the slit as a function of transverse 
position on the slit (ignoring objects outside of the slit, 
as they will in general fall on another slit), with respect 
to an object at the center of the slit, 



T(fi) 



erf [ w+2x 0\ +erf 




2 erf 


so ' _ 





(6) 



Here w is the slit width, s is the seeing full- width at half 
maximum (FWHM) during survey observations, sq is the 
seeing FWHM during standard star observations, x(Sl) 
is the transverse dis tance of the source from the center of 
the slit, and /3 = \/ m ( 2 )- The 

minimum value of T(f2) 
in our survey is about 0.8, so it has a minor effect on the 
computation of Ln m (Cl, z). 

6.2. Volume as a function of source redshift and 
luminosity 

The total volume of our survey sensitive to a source of 
Lyman a line luminosity L is the integral over all volume 
elements in the survey with Lii m (l~2, z) < L, 



V C (L) 



f fdV c {n,z) H[L-L lim (Sl,z)}, 
Jn J z 



(7) 



where H(y) is the step function defined with H(y > 0) = 
1. 

Figure shows the redshift distribution of our survey 
volume as a function of Ln m (Sl, z). The general slight 
decrease toward high redshift is due to the evolution of 
the line element with redshift, and the modulation is due 
to the wavelength-dependent limiting line flux (see Fig. 

EJ. 

We divide our survey arbitrarily into two redshift bins, 
4.6 < z < 5.6, and 5.6 < z < 6.7. There is no natural 
binning choice, but by breaking our survey at z = 5.6 we 
retain almost equal survey volume (at the brightest lumi- 
nosities) in each bin. However this places all three of our 
confirmed high-redshift detections into the 4.6 < z < 5.6 
bin, and none in the 5.6 < z < 6.7 bin. If we had broken 
the bins at z = 5.5, the number density in the lower red- 
shift bin would decrease, and the number density in the 
higher-redshift would increase, that is, the removal or in- 
clusion of a source substantially outweighs the change in 
volume associated with changing the redshift binning. 

In Fig. 1101 we plot the survey volume sensitive to a 
source of luminosity L for each of our two redshift bins 
(represented by the two different symbols). At high lumi- 
nosities there is no dependence of the survey volume on 
luminosity, because sources at such high luminosities are 
so bright that we would detect them at any magnification 
factor or redshift in our survey. At L = 10 42 erg s" 1 the 
high-redshift bin has less volume because of the stronger 
sky lines at longer wavelengths (see Figs. |H] and |!5J) and 
larger luminosity distance compared to the low-redshift 



bin. At lower luminosities the survey volume for both 
bins falls off steadily and similarly. This is a because 
only strongly magnified volume elements contribute to 
the survey volume, and the number of volume elements at 
a given magnification is not sensitive to redshift (Fig.EJ). 

6.3. Number Density 

Our survey detected three confirmed and two likely 
z > 4.5 sources, so to estimate a relatively robust num- 
ber density parameter, and for comparison with other 
surveys, we compute a cumulative number density of 
sources. We construct the cumulative number density 
at each value of the Lyman a line luminosity L (in each 
redshift bin) by evaluating the survey volume at that lu- 
minosity (see above), and then counting the number of 
detected sources brighter than L in the survey volume. 

FigurelTDshows n(> L), the number density of sources 
with Lyman a line luminosities greater than L, for our 
two redshift bins, considering only the three confirmed 
sources. There are only upper limits at L > 10 erg s _1 
because although all three detected sources are in the 
survey volume, none was that luminous. Our most lumi- 
nous source is 7.4 x 10 42 erg s _1 , so the first data point 
appears at L — 10 42 5 erg s _1 (in the low-redshift bin). 
All three of our confirmed detections contribute to the 
L = 10 415 erg s _1 point because all are brighter than 
that limit, and each would have been detected even if 
its luminosity were only 10 415 erg s _1 . In contrast at 
L = 10 41 erg s _1 , all three confirmed detections are still 
brighter than this luminosity, but only one is located in- 
side the L — 10 41 erg s _1 survey volume. At yet fainter 
luminosities we are back to upper limits because none of 
the three confirmed detections would have been discov- 
ered if it were fainter than 10 41 erg s _1 . 

In our high-redshift bin we have no detections, and 
thus can provide only upper limits at all luminosities. It 
is clear that though we can rule out a strong increase in 
the number counts of Lyman a emitters at 5.6 < z < 6.7 
compared to 4.6 < z < 5.6, we cannot further constrain 
the number-count evolution. In particular, our results 
are consistent with no evolution or a decrease with in- 
creasing redshift in the Lyman a source counts as a func- 
tion of redshift over 4.6 < z < 6.7. 

All upper limits and error bars in Fig.^Jare 95% confi- 
dence limits calculated using Poisson statistics. We have 
conceptually divided our survey into sub-surveys sensi- 
tive down to different Lyman a line luminosities, but 
these sub-surveys are not independent (and in fact highly 
correlated). If, for example, a theoretical curve passed 
just through the upper error bars of two points, our data 
would indicate roughly a 95% inconsistency, not a 99.8% 
inconsistency. 

The right and top axes of Fig. ^2 are labeled with unit 
conversions of the left and bottom axes, assuming the 
data fall at z = 5. These serve to allow a rough ref- 
erence of our results to be easily read off in the other 
units commonly used to describe the abundance of Ly- 
man a emission-line galaxies. The right and top axes are 
inapplicable to our high-redshift bin upper limits. 

FigurelTTlshows n(> L), the number density of sources 
with Lyman a line luminosities greater than L, for our 
two redshift bins, considering all five confirmed and likely 
sources. Since we added detections while keeping the sur- 
vey volume fixed, the number densities increased. The 
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number densities of this sample are still marginally con- 
sistent with the 95% confidence limits from Fig. ^2 

7. COMPARISON WITH OTHER OBSERVATIONS AND 
WITH THEORY 

7.1. Comparison with Other Observations 

In Table 01 we present parameters inferred from our 
survey and existing z ~ 5 galaxy surveys. The first 
five entries in the table describe our survey, divided by 
redshift bin and sub-survey limiting Lyman a line lumi- 
nosity. The remaining surveys above the horizontal rule 
are other Lyman a emission-line surveys. Surveys be- 
low the horizontal rule are Lyman-break galaxy surveys, 
described later in this section. 

There are two entries in the number of sources col- 
umn for each row corresponding to our data. The first 
number is the total number of confirmed and likely Ly- 
man a sources in that sub-survey, and the second num- 
ber, in parentheses, in the number of those that are con- 
firmed. We report the corresponding number densities 
analogously in the density column. 

The limiting luminosity, volume, and number density 
of each of the previously published Lyman a emission- 
line surveys appearing in Table^Jdo not always appear in 
the corresponding reference. As necessary we have used 
the published information to calculate those values our- 
selves (for example, converting a limiting line flux and 
redshift into a limiting line luminosity). We expect that 
the final results published by the groups may differ some- 
what. In particular, there may be a publication bias 
toward surveys with discoveries, so it is possible there 
is some bias in the data presented toward higher num- 
ber density. We have included only systematic Lyman a 
emission surveys, because reconstructing the survey vol- 
umes of published serendipitous discoveries was not pos- 
sible. 

FigED plots the data from the 4.6 < z < 5.6 bin of 
our survey (as solid circles) with the data from the other 
Lyman a surveys listed in Tabled (open squares). Note 
that these points are in general from different redshifts, 
and no redshift correction has been applied. The error 
bars shown are 95% confidence limits assuming Poisson 
errors. 

A comparison of our data with published results shows 
that, by utilizing strong lensing, we have provided 
meaningful upper limits on the population of Lyman a 
emission-line galaxies two orders of magnitude fainter 
than previous surveys, in addition to providing confirmed 
data one order of magnitude fainter. With existing ob- 
servational facilities, lensed surveys are the only way to 
probe to such depth. 

At L = 10 42 5 erg s _1 , where our survey overlaps other 
Lyman a surveys, there is marginal consistency between 
our data and published results. Most of the other Ly- 
man a surveys are narrow-band photometric surveys (in 
particular the three surveys with > 10 sources in Ta- 
ble , where the points plotted do not represent con- 
firmed sources, but rather photometric candidates cor- 
rected for the spectroscopic success rate of a small sam- 
ple. As noted previously this figure represents data re- 
ported inhomogeneously, so some of the scatter may be 
related to the different redshift ranges and Lyman a 
equivalent-width criteria of the surveys, as well as er- 
rors in contamination estimation and possible errors in 



our interpretation of published information. 

We detected no sources at z > 5.6. This is marginally 
inconsistent with existing data of source densities at z — 
5.7 and z = 6.5. However, our lack of sources at z > 
5.6 compared to other surveys is qualitatively consistent 
with our smaller number density of 4.6 < z < 5.6 sources 
at L = 10 42 - 5 erg s" 1 . 

For comparison, we have plotted results 
from four z ~ 5 Lyman-break galax y (LBG) 

surveys llStanwav. Bunker, fc McMahonl 2003 ; 

Yan. Windhorst, fe CoherJ l2003t llwata et all bo05 : 



Fontana et all 12003 ILehnert fc Bremerl I2002|) The 



parameters for these surveys are listed in Table 01 below 
the horizontal rule. Again we have converted published 
data into number density as necessary, and taken a 
further step to plot those points on a Lyman a line 
luminosity scale: the LBG survey limit was converted 
into a rest-frame UV continuum limit, th en into a 
star-fo rmation rate using the relation of Kennicutt 



1998), then into a Lyman a line luminosity assuming 
1 Mq yr _1 of star formation prod uces 10 42 erg s^ 1 in 
the Lyman a line l)Kennicuttlll998l after converting Ha 
luminosity into Lyman a luminosity). No unmitigated 
conclusions can be drawn from this comparison, though 
it is intriguing that the z ~ 5 LBG surveys may 
be discovering the same population as the Lyman a 
emission-line galaxies, if the Lyman a line is typically 
1/3 the value expected based on the UV continuum 
SFR. This is similar to the ratio observed in the z = 5.7 
sample of Lyman a emitters of l Aiiki et alJ (12003^ 1 and 
in two galaxie s at z = 6.5 bv IHu et al.1 l)2002b|i and 
iKodaira et alJ (|2003f) . However, four of six confirmed 
4.8 < z < 5.8 galaxies selected by t he Lyman-break 
technique by ILehnert fc Bremerl l|2002j) have Lyman a 
line fluxes less than 10% of the values naively predicted 
from their UV continuum SFRs. 

7.2. Comparison with Theoretical Models 

First we co mpare our results w i th th e Lyman a emit- 
ter model of lHaiman fc Spaansl l)1999fl . who predicted 
the abundance of Lyman a emitters over a range of red- 
shifts and luminosities. In Figs. 1131 and 1141 we plot the 
predictions of their fiducial model at z = 5 as a long- 
dashed curve from log 10 L — 40.5 to 42.5. The shape of 
the lu minosity function predicted bv lHaiman fc Spaansl 
( 1999) is similar to our observed points, but their fidu- 
cial model predicts approximately an order of magnitude 
more sources than we find. Their models could be recon- 
ciled with our data by adopting mass-dependent values of 
the star- format ion efficiency or covering fraction of dusty 
clouds inside the galaxies. 

As a basis for comparing our results with a simple the- 
oretical model, in Figs. 1131 and 1141 we re-plot the lumi- 
nosity function from Fig. ^ (converted into cumulative 
form) , assuming no suppression of the SFR in low-mass 
halos. This simple interpretation of the Lyman a lumi- 
nosity function relates the number density of galaxies to 
dark matter halos. We then converted baryons within 
those dark matter halos into stars, and stellar ionizing 
light i nto Lyman a photons. Unlike lHaiman fc Spaansl 
( 1999) we made no attempt to model the radiative trans- 
fer of the Lyman a photons. 

Specifically, we assumed that 10% of the baryons in 
each halo were converted into stars every halo dynamical 
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time (denned as the ratio of the halo virial radius to the 
halo circular velocity at the virial radius). The Hubble 
time at z — 5 is roughly 10 times longer than the halo 
dynamical time, thus it is possible for such halos to main- 
tain steady star formation at this rate. Star-formation 
rate was converted into an ionizing-phot on rate using a 
Salpe ter IMF with 1/20 solar metallicity (|Leitherer et alJ 
1999). We assumed that 10% of the ionizing photons es- 
cape the emitting galaxy, and that 2 /3 of the remaining 
photons are converted into Lyman a emission. 

The luminosity function predicted by this simple model 
provides a poor fit to our data. In the context of the 
model, it is instructive to think of two modifications 
that would make the predicted luminosity function more 
closely match our data. The first is to decrease the ef- 
ficiency factors used to convert halo mass to Lyman a 
luminosity in a given halo. Alternately, the model curve 
could be brought into agreement with our data if the effi- 
ciency factor for the production of Lyman a was correct 
for a fraction of halos, but the rest had no observable 
Lyman a emission at the time of observation. 

There are three efficiencies that contribute to the over- 
all conversion of halo mass into Lyman a luminosity, 
namely the fraction of baryons converted into stars per 
halo dynamical time, the stellar emissivity of ionizing 
photons, and the fraction of ionizing photons observed 
as Lyman a emission. Lowering the combined efficiency 
by 1.5-2 orders of magnitude would bring the model lu- 
minosity function into close agreement with our data. 
There is some difference in the shape of the curve com- 
pared to our data, but this difference is not significant. 

In the case where only a fraction of halos contain Ly- 
man a emitters, we would require about 1% of halos 
to contain emitters at any given time. This could be 
because, in contrast to the simple model we described, 
star formation is episodic in nature. In addition, there 
could be a timescale associated with the escape of Ly- 
man a photons, such that, for example, dust extinguishes 
Lyman a emission at the beginning of a starburst, but 
eventually the dust is expelled and the Lyman a e mis- 
sion line becomes visible (e.g.. IShaolev e t al. 2003). If 
only some halos contain galaxies, for whatever reason, 
then this formalism of assuming only a fraction of ha- 
los contain Lyman a emitters can also be used, where 
the fraction now represents a filling factor, rather than a 
duty cycle. 

The resolution of the discrepancy between our data 
and the model curve has important implications for the 
mass of the halos that contain the Lyman a emitters. If 
we over-estimated the Lyman a photon production effi- 
ciency in our model, then the association between halo 
mass and Lyman a emitter expressed in Figs. 1X31 and 1141 
is not correct: the halo mass of our population of emit- 
ters at L ~ 10 41 - 5 erg s" 1 should be - 10 11 M . This 
is the largest mass that could be inferred for this pop- 
ulation, assuming a maximum of one Lyman a emitter 
per halo. From the arguments of Section |3 we may ex- 
pect that halos with masses > 10 10 M Q to form stars 
roughly similarly to one another, i.e., though negative 
feedback may be important in regulating star formation, 
it is ineffective in halos this massive. Consequently, for 
this low-efficiency, high halo-mass solution to the discrep- 
ancy, we expect that our data should follow the shape of 
the dark-matter halo mass function, which they do. 



In contrast, if we resolve the discrepancy between the 
model and our data by assuming that the efficiency we 
calculated is correct for a fraction of halos, and the rest 
are empty of Lyman a emission, then the mass associ- 
ation in Figs. 1131 and 1141 is correct. This implies a halo 
mass of only ~ 10 9 ' 5 M a for our Lyman a emitters at 
L ~ 10 415 erg s~ x . Depending on the characteristic mass 
scale where negative feedback becomes a dominant pro- 
cess, the Lyman a luminosity function may already de- 
viate from the shape of the mass function at ~ 10 95 M Q 
(see Fig.0. Our data are slightly flatter than the mass 
function, and consistent with any of the luminosity func- 
tions plotted in Fig. ^ 

Our theoretical interpretation so far has relied exclu- 
sively on our data, which is consistent with the shape of 
the relevant halo mass function. However, if we consider 
all of the available Lyman a data, there is some evidence 
for a flatter Lyman a luminosity function. The hetero- 
geneous nature of the Lyman a survey data plotted in 
Figs. El an d El makes it difficult to draw firm conclu- 
sions, but a combination of data at L > 10 42 5 erg s _1 
suggests that our data point at L — 10 415 erg s" 1 may 
be 0.5-1 dex lower than an extrapolation of the Lyman a 
luminosity function from higher luminosity, assuming the 
luminosity function shape matches the mass function 
shape. Thus we conclude that our data, in combina- 
tion with other Lyman a surveys, suggest that strong 
negative feedback is suppressing the star-formation rate, 
a nd thus Lyman a lum inosity, in our sources. 

IHamana et alJ 1120 03) used clustering data to estimate 
the mass of the halos containing Lyman a emitters at 
z = 4.9. They concluded that the c haracteristic halo 
mass of those sources is 5 x 10 12 M Q l|Shimasaku et alJ 
120031 find a halo mass of - 10 12 M Q for similar z = 4.9 
emitters on the basis of a large-scale structure feature in 
their survey). This conclusion would support the low- 
efficiency, high halo-mass solution to the difference be- 
tween our model luminosity function and our data. How- 
ever, the number density of z = 4.9 emitters is larger, by 
about a factor of five to ten, than the number density 
of 10 12 M halos (jHamana et al.M2003|) . This implies, 
contrary to our assumption above, that there is more 
than one Lyman a source per halo. The virial radius of 
a z — 5, 10 12 Mq halo is 8.5", so multiple sources in- 
side a single halo may be observed as separate sources, 
though this should create a very distinct signature in the 
spatial distribution of sources (or extended nature, if the 
sources are unresolved) that has not been reported by 
other Lyman a emitter surveys. 

While current information on the masses of Lyman a 
emitter halos is still limited, progress will continue to be 
made at L > 10 42 5 erg s _1 by large Lyman a surveys. 
Unfortunately, surveys for low-luminosity Lyman a emit- 
ters will not provide sufficient survey area for clustering 
studies in the near future. Lensed surveys such as ours, 
in particular, do not lend themselves to easy clustering 
analysis, because the contiguous survey volume is very 
complex and limited in size by the mass of the lensing 
foreground cluster. As an aside we comment that Poisson 
errors dominate the uncertainty in the number densities 
plotted in Figs. El and El (and Figs. El and El for our 
survey), assuming the maximum mass for the halos con- 
taining our emitters (i.e., every halo contains a source; 
see above for a caveat), and using the clustering formal- 
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ism of IMo & White) 1I2OOI . 

Until the advent of large-area, low-luminosity Lyman a 
surveys, the only constraint on the mass of the halos 
containing the emitters, and thus the only path toward 
understanding the suppression of star-formation in low- 
mass halos, lies in detecting source populations with high 
number densities, such that the halo mass function, and 
the assumption that there is at most one source per halo, 
can be used to infer a maximum halo mass for the popu- 
lation of Lyman a emitters. This is strong motivation for 
future surveys to continue to use strong lcnsing to survey 
small volumes to considerable depths for faint, Lyman a 
emitting sources. 

8. SUMMARY 

We performed a systematic survey for Lyman a emis- 
sion at 2.2 < z < 6.7 using strong lensing from 
intermediate-redshift clusters of galaxies to boost our 
survey sensitivity to unprecedented depths. We detected 
three confirmed and two likely Lyman a emitting galax- 
ies at 4.7 < z < 5.6, with Lyman a line luminosities of 
2.8 x 10 41 erg s _1 < L < 7.4 x 10 42 erg s -1 . Our survey 
covered 4.2 arcmin 2 on the sky, with a maximum vol- 
ume of 4 x 10 4 Mpc 3 over 4.6 < z < 6.7. We find no 
evidence for redshift evolution of the number density of 
Lyman a emitting galaxies between z ~ 5 and z ~ 6, 
though our data are also consistent with a decrease in 
number density with increasing redshift. 

We present the first meaningful constraints on the the 
luminosity function of Lyman a emitters at 4.6 < z < 5.6 
over the Lyman a luminosity range 10 40 erg s _1 < L < 
10 42 erg s _1 , corresponding to inferred star- formation 



rates of 0.01-1 M© yr _1 . From a consideration of the 
number density of dark-matter halos, we conclude that 
our population of sources at L ~ 10 415 erg s^ 1 resides 
in halos of mass < 10 11 M©. 

Our number density data are consistent with a Ly- 
man a luminosity function with the same shape as the 
halo mass function, but a consideration of all available 
Lyman a survey data implies that we have observed a 
flattening of the Lyman a luminosity function with re- 
spect to the halo mass function. We may have detected 
evidence of the suppression of star-formation in low-mass 
halos at high redshift, as predicted by theoretical models 
of galaxy formation. 
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Fig. 1. — The star-formation rate function based on a dark matter halo model. The solid curve is the abundance (on the right axis) of 
halos at z = 5 as a function of the mass on the top axis. The bottom axis is a simple conversion of the halo mass into the expected star 
formation in that halo, so that the solid curve may also be read as a theoretical star-formation rate function using the bottom and left 
axes. The broken curves represent cases where star formation in lo w-ma ss halos has been suppressed; each is labeled with the logarithm of 
the mass scale below which suppression occurs. See Sections 151 and IY. 21 for details. 



Table 1 
Clusters surveyed 



Cluster 


Redshift RA 


A 


Dec b 


Lens Model Reference 


Abell 68 


0.255 


00 


36 


59 


+09 09 


(1) 


Abell 370 


0.375 


02 


37 


18 


-01 48 


(2) 


Abell 773 


0.217 


09 


14 


30 


+51 55 


(1) 


Abell 963 


0.206 


10 


17 


09 


+39 01 


(1) 


Abell 1689 


0.183 


13 


09 


00 


-01 06 


(3) 


C11358. 1+62.45 


0.328 


13 


59 


51 


+62 31 


(4) 


Abell 2218 


0.176 


16 


35 


42 


+66 19 


(5,6) 


Abell 2219 


0.226 


16 


38 


54 


+46 47 


(1) 


Abell 2390 


0.228 


21 


53 


35 


+17 40 


(7) 



a units of HH MM SS 
b units of +DD MM 



Table 2 
LRIS Survey Observations 



Date 


Cluster 


Position Angle 3, 


Integration timc b 


Photometric? 


Mar 2000 


Abell 773 


-46.8 


20 


Yes 




Abell 1689 #1 


84.1 


23 


Yes 


Apr 2001 


Abell 1689 #2 


43 


10 


Yes 




Abell 2218 #1 


-44 


10 


Yes 




C11358 


-15 


12 


Yes 


Oct 2001 


Abell 370 #1 


-8 


14 


Yes 


Apr 2002 


Abell 963 


3.6 


14 


No 




Abell 2218 #2 


-49.2 


14 


No 


May 2002 


Abell 1689 #3 


12.3 


20 


No 




Abell 2219 


-69 


14 


No 


Sep 2002 


Abell 370 #2 


1.7 


14 


No 




Abell 2390 


-63 


12 


No 




Abell 68 


-40 


12 


No 



a in degrees North through East 
b in ksec 
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Fig. 2. — Survey geometry for a selected sample of our clusters. For each cluster, superposed on the HST /WFPC2 image are the critical 
lines for a source at z = 5 (dotted lines). The solid curves bound areas where the magnification for such a source exceeds a factor of 10. 
The regions bounded by parallel straight lines are the long-slit survey area. Numeric labels correspond to the key in Table |21 The axes are 
labeled in arcseconds. (upper left) Abell 68. (upper right) Abell 370. (middle left) Abell 773. (middle right) Abell 963. (lower left) Abell 
1689. (lower right) Abell 2218. 
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Fig. 0cont.— (upper left) Abell 2219. (upper right) Abell 2390. (lower left) C11358. 



Table 3 

Lyman a Emission-Line Candidates 



Cluster ID 


RA 


Dec 


A a 


z h 


Flux c 


Comments 




Abell 370.1.f 


2:39:50.60 


-1:33:45.0 


4628 


2.80 




confirmed (LRIS] 




Abell 370. l.g 


2:39:51.80 


-1:35:57.6 


4630 


2.80 




confirmed (Ivison et al. 


2002) 


Abell 963.1.cd 


10:17:05.10 


39:03:30.5 


5191 


3.27 




confirmed (ESI) 




Abell 963.1. efg 


10:17:04.77 


39:03:11.0 


5617 


3.62 




confirmed (ESI) 




Abell 1689. 2.f 


13:11:25.38 


-1:20:52.4 


7141 


4.82 


3.0 


confirmed (Frve et al. 


2002) 


Abell 2218. 1.a2 


16:35:51.75 


66:12:45.6 


8001 


5.58 


4.4 


confirmed ("ESI. Ellis et al. 20011 




16:35:51.89 


66:12:51.5 








2nd image 




C11358.1.ef 


3:59:49.19 


62:30:44.8 


7205 


4.92 


10 


confirmed (Franx et al. 




Abell 773. l.c 


9:17:55.31 


51:44:26.6 


6978 


4.74 


1.1 


likely 


Abell 963. l.d 


10:17:04.45 


39:01:47.1 


7025 


4.77 


0.69 


likely 




Abell 2218. l.al 


16:35:45.25 


66:13:26.4 


4216 


2.47 




likely 




Abell 2218.2. b 


16:35:48.78 


66:12:24.9 


3928 


2.23 




likely 





a Wavelength of emission line in units of A 

b Source redshift assuming emission line is Lyman a 

c Observed line flux (uncorrected for lensing, corrected for transparency) in units of 10 — 17 erg s — 1 cm -2 
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Fig. 3. — Two-dimensional spectra of the three confirmed z > 4.5 galaxies detected in the survey. The spectra, from left to right, are the 
z = 4.89 source in Abell 1689, the z = 5.57 source in Abell 2218, and the z = 4.92 source in C11358. The wavelength coverage in all the 
spectra is 6800 to 8430 A, increasing bottom to top. 
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Fig. 4. — Two-dimensional spectra (left) and HST images (right) for our three confirmed Lyman a sources. In the spectra the horizontal 
axes are labeled in A, and the vertical axes in arcsec; the image axes are labelled in arcsec. (upper) Abell 1689. 2. f. (middle) Abell 2218. 1.a2. 
(lower) C11358.1.ef. 




. r h 




Fig. 5. — Two-dimensional spectra (left) and HST images (right) for our two likely Lyman a sources. In the spectra the horizontal axes 
are labeled in A, and the vertical axes in arcsec; the image axes are labelled in arcsec. (upper) Abell 773. I.e. (lower) Abell 963. l.d. 
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Fig. 6. — Magnification due to the cluster Abell 2218 of background sources, as a function of position and redshift. This figure shows 
the magnification factor at positions along a 175" longslit at one of our survey positions in Abell 2218 (within pointing "#1" in Fig. 
The solid curve is for sources at z = 4.3, and the dotted curve is for sources at z = 6.8. 




100 1000 
.iiagnif ication 



10000 



Fig. 7. — Cumulative histogram of the magnification factor over the entire survey area. The solid and dotted curves show the magnification 
histograms for sources at z = 4.3 and z = 6.8, respectively. 
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Lyman a liedshiFl 



5 5.5 6 6.5 




7000 7500 8000 8500 9000 
Wavelength [A] 

Fig. 8. — Limiting Lyman a flux as a function of wavelength for the survey exposure time of 2000 s, assuming photometric conditions. 
The curve is the 5-cr flux limit to detect an emission line in a 1.3" by 7.7 A aperture at that wavelength. The top axis is labeled with the 
redshift corresponding to observed Lyman a falling at the wavelength on the bottom axis. The strong fluctuations in /ii m (z) are caused 
by atmospheric emission lines, and the rise at high redshift is due to decreased instrumental sensitivity. 




Redshift (z) 

Fig. 9. — The redshift distribution of our survey volume. The curves, from top to bottom, represent the redshift distribution of subsurveys 
down to limiting luminosity log 10 L = (43.5, 43, 42.5 , 42, 41.5); for yet lower values of the limiting luminosity, the curves have the shape of 
the bottom curve, but are scaled down (see Fig. 1101 . 
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Fig. 10. — Survey volume sensitive to sources with Lyman a line luminosity L. The points show the volume of the survey within which 
a source of luminosity L would be detected. The survey has been divided into two redshift ranges: the volume with 4.6 < z < 5.6 is shown 
with squares, and the volume with 5.6 < z < 6.7 is shown with diamonds. 
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Fig. 11. — Number density of sources brighter than luminosity L, for the three confirmed z > 4.5 sources. Each luminosity is treated as 
an independent sub-survey down to luminosity limit L, and the cumulative number density is calculated from the number of sources in the 
sub-survey and the volume of the sub-survey. The survey was divided into two redshift bins: the solid lines and points are for 4.6 < z < 5.6, 
and the dotted lines (off-set slightly for clarity) are for 5.6 < z < 6.7. The error bars are 95% limits calculated using Poisson statistics. 
The top and right axis are labeled in units for comparison of the 4.6 < z < 5.6 bin results with other work: the left and bottom axes were 
transformed assuming a redshift of z = 5.0. 



Low-luminosity Lyman a sources at high redshift 
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Fig. 12. — Number density of sources brighter than luminosity L, for the three confirmed plus two likely z > 4.5 sources. As in Fig. 1111 
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Fig. 13. — Number density of z ~ 5 sources brighter than L, for Lyman q surveys and 2 ~ 5 Lyman-break galaxy (LBG) surveys. 
The solid circles are our cumulative number densities of confirmed sources brighter than L for sub-surveys within the 4.6 < z < 5.6 bin 
(described in Fig. 1111 . The open squares are the cumulative number densities of sources brighter than L inferred from other z ~ 5 Lyman a 
line surveys, and the crosses are data from z ~ 5 LBG surve ys. The LBG surveys were converted to equivalent Lyman a line luminosities 
(see Section 17. II . The long-dashed curve is a prediction from Haiman & SDaans] I199HI) . The solid curve is the cumulative number density 
of halos above the mass given on the top axis; the vertical scale is the same. The data are described in Table 151 
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Fig. 14. — Number density of z ~ 5 sources brighter than L, for Lyman a surveys and z ~ 5 Lyman-break galaxy (LBG) surveys. The 
solid circles are our c umu lative number densities of confirmed and lik ely sources brighter than L for sub-surveys within the 4.6 < z < 5.6 
bin (described in Fig. EU- The other symbols are the same as in Fig. 1131 



Table 4 
Galaxy surveys at z ~ 5 



Redshift 


logic, L a 


Number 13 


Volumc c 


Dcnsity d 




Reference 


4.6-5.6 


40.5 


1(0) 


0.0053 


1.8 x 10~ 2 (0) 




(1) 


4.6-5.6 


41 


4(1) 


0.037 


1.1 x 10~ 2 (2.7 x 10" 


" 3 ) 


(1) 


4.6-5.6 


41.5 


3(3) 


0.20 


1.5 x 10~ 3 (1.5 x 10" 


" 3 ) 


(1) 


4.6-5.6 


42 


3(2) 


0.92 


3.3 x 10~ 2 (2.2 x 10" 


4 ) 


(1) 


4.6-5.6 


42.5 


1(1) 


2.0 


5.1 x 10" 4 (5.1 x 10" 


-4) 


(1) 


5-6 


42.28 


1 


0.018 


5.6 x 10~ 3 




(2) 


5.7 


42.61 


1 


0.55 


1.8 x 10~ 4 




(3) 


6.5 


42.95 





6.1 







(1) 


6.5 


42.35 


1 


0.012 


8.5 x 10" 3 




(4) 


1.9 


42.46 


87 


9.2 


9.5 x 10~ 4 




(5) 


5.7 


42.81 


13 


28 


4.6 x 10~ 5 




(6) 


6.5 


42.89 


16 


20 


8.0 x 10~ 6 




(7) 


5.7 


42.95 


6 


13 


4.6 x lO" 5 




(8) 


5.7 


42.72 





0.11 







(9) 


5.7 


42.40 





0.073 







(9) 


~ 5.8 


43.35 


6 


18 


3.3 x 10~ 5 




(10) 


~ 6.3 


43.32 


26 


1.7 


1.5 x 10" 3 




(11) 


~ 5 


42.96 


10 


1.4 


7.1 x 10" 4 




(12) 


~ 5.5 


43.40 


2 


2.4 


8.3 x 10~ 5 




(13) 


~ 5.5 


43.28 


16 


7.0 


2.3 x 10~ 3 




(13) 


~ 5.3 


43.18 


6 


12 


5.0 x 10~ 5 




(14) 



a Survey limiting L in units of erg s — 1 
b Number of sources detected in survey 
c Volume of survey in units of 10 4 Mpc 3 
d Number density of sources, in units of Mpc -3 



